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In the next decade gravitational waves might be detected using a pulsar timing array. In an effort to de- 
velop optimal detection strategies for stochastic backgrounds of gravitational waves in generic metric theories 
of gravity, we investigate the overlap reduction functions for these theories and discuss their features. We 
show that the sensitivity to non-transverse gravitational waves is greater than the sensitivity to transverse grav- 
itational waves and discuss the physical origin of this effect. We calculate the overlap reduction functions 
for the current NANOGrav Pulsar Timing Array (PTA) and show that the sensitivity to the vector and scalar- 
longitudinal modes can increase dramatically for pulsar pairs with small angular separations. For example, the 
J 1853+ 1303-1 1857+0943 pulsar pair, with an angular separation of about 3°, is about 10'' times more sensitive 
to the longitudinal component of the stochastic background, if it's present, than the transverse components. 



I. INTRODUCTION 

General relativity is among the most successful theories of 
physics in the 20th century, passing all current weak-field, 
slow motion tests with flying colors. Progress in cosmology 
and high energy physics over the course of the last fifty years, 
however, has brought with it questions that may be unanswer- 
able in the context of general relativity. The accelerated ex- 
pansion of the universe, the dark matter problem, and infla- 
tion have led some authors to re-examine general relativity 
and attempt to modify it to explain some of these puzzles. Ad- 
ditionally, the incompatibility between general relativity and 
quantum field theory may be an indication that modifications 
to general relativity are necessary. 

A number of alternative theories of gravity have been pro- 
posed to address some of these problems. Those which satisfy 
the Einstein Equivalence Principle are called metric theories 
of gravity. In these theories, the only gravitational fields that 
may influence matter are the components of the metric tensor 
(jfpix. Additional fields play the role of generating spacetime 
curvature. Metric theories are grouped broadly into several 
categories: scalar tensor theories, in which a dynamical scalar 
field is present in addition to the metric (see Refs. lHJIH); 
vector-tensor theories, which contain a dynamic gravitational 
four-vector field in addition to the metric (see Refs. |7, 9 - 12 1); 
and bimetric theories, which are characterized by "prior" ge- 
ometry contained in dynamical scalar, vector or tensor fields 
(see Refs. 10 IS [13). 

Gravitational wave astronomy promises not only to open a 
new observational window on the universe, but also to pro- 
vide a new testing ground for general relativity. In a general 
metric theory of gravity, the six independent components of 
the Riemann tensor provide up to six possible gravitational 
wave (GW) polarization states, four more than those allowed 
in general relativity. Detection of any extra GW polarization 
states would be fatal for general relativity. A non-detection 
could be used put constraints on the parameters of alternative 
theories of gravity. 

Several international efforts are currently underway to de- 
tect GWs. Of these the most promising on the 5-10 year 
timescale are ground-based laser interferometers 1 14] and pul- 



sar timing arrays 1 15 1, which aim to detect GWs in the lO-lO"^ 
Hz and 10^^-10^^ Hz ranges, respectively. Potential sources 
for low frequency GWs (10^^-10^^ Hz) include binary su- 
permassive black hole mergers il6l , cosmic superstrings 1 17|, 
relic gravitational waves from inflation ifTSl . and a first order 
phase transition at the QCD scale lfT9]| . 

Previous work on stochastic backgrounds of gravitational 
waves in the context of alternative theories of gravity has 
shown that three ground-based interferometers are sufficient 
to disentangle the polarization content of a general metric the- 
ory of gravity |20|. For pulsar timing arrays the form of the 
correlation between pulsar pairs as a function of pulsar pair 
angular separation depends on the polarization content of the 
theory [21 1. Additionally it has been shown that pulsar tim- 
ing arrays have a greater sensitivity to longitudinal and vector 
polarization modes than to transverse modes |2T] Ell- 
in this paper we investigate the problem of stochastic GW 
detection using PTAs in the context of the optimal statistic. 
We compute the expected cross correlations for pulsar tim- 
ing arrays for the case of stochastic backgrounds of GWs for 
any metric theory of gravity. The expected cross correlations 
are proportional to the so-called overlap reduction function, a 
geometrical factor that captures the loss of sensitivity due to 
detectors not being co-located or aligned. We explain various 
features of the overlap reduction functions including the phys- 
ical origin of the increased sensitivity to scalar-longitudinal 
and vector polarization modes. In Section [H] we use a coor- 
dinate independent approach to describe the redshift of pulsar 
signals from passing GWs. In Section III we write the opti- 



mal cross-correlation filter by maximizing the signal to noise 
for a pulsar pair, and define the overlap reduction function for 
GWs of any metric theory of gravity. In Section|lV]we discuss 
the effect of GWs of various polarizations on the pulsar-Earth 
system, and the physical origin of the increased sensitivity to 
longitudinal and shear modes. This effect is most easily un- 
derstood in the frequency domain. In Section |Vj we write 
down explicitly the form of the overlap reduction function for 
transverse GWs and discuss the form of the function for non- 
transverse GWs. We find that for the scalar-longitudinal and 
vector (shear) modes, the overlap reduction functions are fre- 
quency dependent in the ranges of frequencies and distances 
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FIG. 1 : Pulsar positions are given with respect to the Solar 
System barycenter (located at the origin). Here 6 and (/> are 
the typical polar and azimuthal angles (as projected from the 
position of pulsar 1), and pulsar 1 and pulsar 2 are separated 
by angle ^. A gravitational wave, characterized by 
polarization angle propagates along the Cl direction. 



relevant to pulsar timing. This is not the case for the transverse 
tensor and breathing modes. 



In Section VI we compute the 
overlap reduction functions for the current NANOGrav PTA 
and show that sensitivity to the scalar-longitudinal and vector 
(shear) modes increases by at least an order of magnitude for 
nearby pulsar pairs for vector modes, and about four orders 
of magnitude for the longitudinal mode. We summarize our 



results in Section Vll Throughout we work in units where the 
speed of light c = 1. 



II. DETECTING GRAVITATIONAL WAVES WITH A 
PULSAR TIMING ARRAY 

The radio pulses from pulsars arrive at our radio telescopes 
at very steady rates. Pulsar timing experiments exploit this 
regularity. Fluctuations in the time of arrival of radio pulses, 
after all known effects have been accounted for, might be due 
to the presence of a GW background. If a GW is present the 
signal from the pulsar can be red-shifted (or blue-shifted). For 
a GW propagating in the direction Q,, the redshift of signals 
from a pulsar in the direction p is given by | 23, 24l 



2[i + n-p] 



(1) 



where /i^ is the metric perturbation and tp, represent the 
times the pulse was emitted at the pulsar and the time received 
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FIG. 2: Motion of test masses in response to GWs of the six 

polarization modes. The plus (+), cross (x), and 
scalar-breathing (6) mode GWs are transverse, while the two 

vector modes {x, y) and the scalar-longitudinal (0 mode 
GWs are non-transverse. Figure reproduced from Nishizawa 
et al. |i20| with permission. 



at the Solar System barycenter, and we have defined 



(2) 



Note that this is the opposite of the sign convention normally 
used in the literature 1231 . Modified gravity theories extend 
the possible polarization modes of GWs present in general rel- 
ativity - the plus (+) and cross ( x ) modes- to a maximum of 
six possible modes. For the two pulsar-Earth system shown 
in Fig. [T] the GW coordinate system is given by 



Cl = (sin 9 cos <j), sin 9 sin (p, cos 9) 
rh — (sin 0, — cos <f), 0) 
n — (cos 9 cos (/), cos sin 0, — sin 9) 



(3) 



where, relative to ||20ll . we have fixed the GW polarization 
angle -if) = — 7r/2 to agree with the conventions in [25J. From 
(|3]l, the GW polarization tensors can be constructed 1261 



e,,-'- = m (g) TO — n (X) e, „■ = m (g) ri + n (g) to 
''J 



''J ' U 



(4) 



where g) is the tensor product and Cl is the direction of GW 
propagation. Here, x and y correspond to the vector (spin- 
1) polarization modes while b and I correspond to the scalar 
(spin-0) breathing and longitudinal modes, respectively. The 
plus, cross and breathing modes are characterized by trans- 
verse GW propagation, while the longitudinal and vector (or 
shear) modes are non-transverse in nature (see Fig.|2]i. 
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Defining the antenna patterns as 



2ii + n-p) 

the Fourier transform of ([T]i may be written as 



(5) 



-l)j2hAif,n)F^{Cl) (6) 



where the sum is over all possible GW polarizations: A = 
+,x,x, y, b, I, and L is the distance to the pulsar The actual 
quantity measured in pulsar timing experiments is the timing 
residual, which is defined as the difference between the actual 
and expected time of arrival (TOA) of a pulse: 



R{t) ^ TOA 



actual 



TOA expected- 



(7) 



The expected TOA for a pulse is modeled and includes daily 
and yearly motion of the Earth, the position and proper motion 
of the pulsar, motion about a binary companion (if applicable), 
etc. The timing residual can be obtained by integrating the 
redshift in time I 



III. GW DETECTION STATISTIC 

In this section we introduce the optimal cross correlation 
statistic ll24l|25 | for stochastic background searches. The op- 
timal cross-correlation statistic involves the calculation of the 
overlap reduction function, a geometrical factor that charac- 
terizes the loss of sensitivity due to detectors not being co- 
located or aligned. We will show how the overlap reduction 
function is computed for non-transverse modes. We follow the 
analysis (for General Relativity) of Allen and Romano f25 1. 

A stochastic background of GWs is produced by a large 
number of weak, independent, unresolvable sources. The en- 
ergy density of this background per unit logarithmic frequency 
is given by 



^gw (/) 



1 



Pcritical din f 



(8) 



where dpg^ is the energy density of the gravitational waves 
and /Ociiticai is the critical energy density required to close the 
universe. 



/^critical 



3Hl 
SttG 



(9) 



where Hq is the Hubble constant. 

The characteristic strain spectrum, hdf), is typically given 
a power-law dependence on frequency so that 



hcif) 



A 



I 



yr- 



(10) 



It may also be expressed in terms of the energy density of the 
background per unit logarithmic frequency, f7gw(|/|): 



hlif) 



2^2 J2"gw 



(l/D- 



(11) 



For an isotropic stochastic background of GWs, the signal 
appears in the data as correlated noise between measurements 
from different pulsars. The i*'' data set is of the form 



,{t) = Zi{t)+n,{t) 



(12) 



where Zi {t) corresponds to the unknown GW signal and ni {t) 
to noise (assumed in this case to be stationary and Gaussian). 
We will make the following assumptions about the expecta- 
tion values of the ni(t): 

{n,{t)) = (13) 
(n,(t)n,(i)) = (14) 

for data sets 

The cross-correlation statistic is defined as 

\ dt [ dt's,{t)sj{t')Q{t-t') (15) 

J-T/2 J-T/2 

where Q{t — t') is the filter function. The optimal filter is 
determined by maximizing the expected signal-to-noise ratio 



SNR = 



(16) 



Here /i is the mean (S) and a is the square root of the variance 
In the frequency domain ^T5\ becomes 



OO pOO 



df dfST{f-ns*{f)SjinQ{n, (17) 

J — OO J —OO 

and the mean fi is 

/OO /* OO 

df / df 6T{f ^ n{s*{m{r))Q{fm 
-OO J —OO 

where 5t is the finite time approximation to the delta function 

sin 7r/i 



<5t(/) 



Assuming that the background is unpolarized and isotropic, 
and stationary. 



(19) 



X sif ~ n\f\-^ngU\f\) 



we obtain 



R*(/)%(/')> 



3Hi 1 



-^s{f - r)\f\-' (20) 



327r3 /3 

X i]gw(|/|)r(|/|) 



Here /? is a normalization factor and we define f24\ 

r(|/|) = I dnie'^'^'f^-'-^+^'^-P'^ - I) (21) 
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FIG. 3: Antenna patterns for plus/cross (a), breathing (b), vector- x/vector-y (c), and longitudinal (d) polarization modes, as 
computed from (|5|. Note that the cross and vector-y modes are identical to plus and vector-x, respectively, but rotated by 45 
degrees and thus do not appear separately here. In this figure, the GW propagates in the positive z-direction with the Earth at 
the origin, and the antenna pattern depends on the pulsar's direction, specified by polar angle 6p and azimuthal angle (j)p. Exact 
expressions corresponding to each figure may be found in Appendix [A] ( |A20| i for the plus mode, ( |A26| l for the breathing mode, 

( |A24| i for the vector-x mode, and ( |A12| i for the longitudinal mode. Note that fixing the GW propagation direction while 
allowing the pulsar position to change is analogous to fixing the pulsar position while allowing the direction of GW propagation 
to change (there is an inherent degeneracy in the GW polarization angle and the pulsar's azimuthal angle (j>p). 



where the sum is over all possible GW polarizations, and the The optimal filter is given by Il24ll25l . 
exponential phase terms correspond to the pulsar term in the q / n 

..n»d„n,..„. «<^>«Mi)' 

where Pi{f) and Pj{f) are the power spectra for the ith and 
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jth pulsar redshift time series that are being cross-correlated 
(see Eq.[T5]l. 

In general relativity, for the frequency and distance ranges 
appropriate to pulsar timing experiments (i.e. for / ^ 1/i), 
the overlap reduction function r(/) approaches a constant 
which is only a function of the angular separation between 
the two pulsars. This constant is proportional to the value 
of the Hellings-Downs curve for the angle between the pul- 
sars ll24l l27l . We will see that for longitudinal modes and 
for tensor modes the overlap reduction function remains fre- 
quency dependent, even for / ^ 1/L, and is considerably 
larger than for the transverse modes. This indicates an in- 
creased sensitivity to such modes. To understand the physical 
origin of the increased sensitivity we first discuss the effect of 
GWs in the more simple case of a single pulsar-Earth baseline. 



IV. GW INDUCED REDSHIFT ON THE PULSAR-EARTH 
SYSTEM 

In this section we will study the redshifts induced by GWs 
of different polarizations on the pulsar-Earth system. From 
(|6]l, the redshift induced by this GW may be written as 



-27ri/L(l+r2-p) 



2{i + n-p) 



(23) 



The factor of 1/2(1 + f2 • p) comes from the relationship be- 
tween the affine parameter A and time t (see Eq. ( |A9| l), and 

hA = hA{f,ti). 

In the region where the GW direction, Ct and the pulsar 
direction, p are anti-parallel, ( |23] l appears to become singular 
due to the 1 + fl ■ p term in the denominator (note that the 
derivative of Ha with respect to the affine parameter vanishes 
in this limit; see ( |A9| l). There is in fact no divergence in the 
redshift induced. In this regime the exponential can be Taylor 
expanded and the 1 + Cl ■ p term in the denominator cancels. 

A Taylor expansion of ( |23| ) can be performed in two cases. 
In the first, when fL <C 1, the metric perturbation is the same 
at the pulsar and at the Earth. This case is often referred to as 
the long wavelength limit. In the second, when 



1 + n-p-^ 



1 



the pulse's direction of propagation and the GW are nearly 
parallel (i.e. the GW is coming from a direction near the pul- 
sar). In this case the metric perturbation at the pulsar when the 
pulse is emitted, and on Earth when the pulse is received, are 
also nearly the same. This is often described in the literature 
in terms of the pulse "surfing" the gravitational wave. 

The surfing description, combined with Eq. ([T}, might 
lead one to incorrectly conclude that the effect of the GW 
should cancel in this case because the metric perturbations at 
the Earth and the pulsar are the same, despite the divergent 
+ Q, ■ p) term in the redshift. In fact, a delicate cancella- 
tion occurs with the divergent term in the denominator which 



is only manifest in the frequency domain. Let the pulse direc- 
tion and the gravitational wave direction be nearly parallel so 
that Vl ■ p — — 1 + (5, where (5^1. Then 



'TTifLp'p'efjhA- 



(24) 



The redshift is proportional to fL, but for finite S increases 
only to the point where the argument of the exponential in 
( |23| ) can no longer be Taylor expanded, at which point it be- 
comes an oscillatory function of fL. Whether the redshift 
is finite in the (5 — > limit depends on the projection term 
p'^p'efjhA- As we will see, the vanishing contribution for the 
tensor modes of general relativity occurs solely because of the 
transverse nature of these waves, and is unrelated to the "surf- 
ing" effect. For longitudinal modes the projection term does 
not vanish, and the increase in sensitivity to such modes orig- 
inates from GWs that come from directions near the pulsar 
To better understand this, we will look at the behavior of the 
redshifts induced by GWs of various modes. 

The redshift for a longitudinal mode GW perturbation is 



cos^ 9 



2(1 + cos ( 



-27ri/L(l+cos9) 



l)hi, (25) 



while the redshift for a plus mode GW perturbation is 



sm 



2(1 + cos ( 



-27rifL{l+cos0) 



l)h+.(26) 



Here we note that the geometrical factor in the redshift for the 
transverse breathing mode differs from ( |26| l only by a sign, 
and our analysis of (|26]l applies equally to the breathing mode. 
In Fig. 4]we plot the geometrical and phase factor \ z{f, Cl) /h\ 
for both the +-mode and the longitudinal mode. We plot these 
for a value of fL in the long wavelength limit (fL = 10~^), 
and for a value in the regime of pulsar timing experiments 
(fL = 10). In the regime of pulsar timing experiments the 
sensitivity is largest for GW directions near the pulsar 6* « tt 
for both polarizations. Although we do not show it here the 
same is true for all other polarization modes. In the longth 
wavelength limit, fL ^ 1, the pulsar-Earth system is most 
sensitive to +-mode GWs coming from the equator, and lon- 
gitudinal GWs from the poles. 

As discussed above, these redshifts appear to become sin- 
gular when 9 —i' TT, but the pulsar term may be Taylor ex- 
panded. Let 9 = TT — S, where (5^1. Then 



Siif,n)^nifLil~S^)hi 



for the longitudinal case, while 



S+{f,n)^7:ifL6^h+ 



(27) 



(28) 



for the plus mode. In the limit as 5 — > 0, z+ vanishes while 
zi becomes proportional to fL. The vanishing redshift of 
is therefore due to the transverse nature of the mode, and does 
not occur for zi, even though in both cases the pulse is "surf- 
ing" the GW. In the time domain, in the « tt region, the 
redshift for both modes goes as 



zi^+{t, Cl) cx Lhi, 



(29) 
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FIG. 4: (color online) Plots of |z(/, d)/h\ for the +-mode (dashed blue) and the longitudinal mode (solid red). We show these 

for JL — 10"^ (a), a value of fL in the long wavelength limit, and (b) JL = 10, a value of fL typical of pulsar timing 
experiments. In the regime of pulsar timing experiments the sensitivity is largest for GW directions near the pulsar 6 -k for 
both polarizations. In the longth wavelength limit, fL <C 1, the pulsar-Earth system is most sensitive to +-mode GWs coming 

from the equator, and longitudinal GWs from the poles. 



One may readily identify the right hand side of p9| l as a ve- 
locity. The interpretation of this result is that, in this limit, the 
redshift is proportional to the relative velocity of the pulsar- 
Earth system. The velocity of the pulsar when the pulse is 
emitted in this limit is approximately equal and opposite to 
the velocity of the Earth when the pulse is received. 

An identical analysis for the shear GW modes produces 
analogous results. Starting from (|6|, the redshift for the 
vector-y mode goes as 



cos 9 sin 



-2TTifL{l+cos( 



(1 + cosB 
The small S expansion yields 

2y{f,n) ^ -2Tn fLS ( 1 



l)hy.(30) 



(31) 



To summarize, one can see that the surfing effect does 
not lead to a vanishing response of the pulsar-Earth system 
to GW waves coming from 9 = n. For the tensor and 
scalar-breathing modes, it is the transverse nature of GWs 
that is responsible for the vanishing response. For the scalar- 
longitudinal modes the response does not vanish — in fact, the 
response increases with both frequency and pulsar distance. 
For the vector modes the response does vanish, but more 
slowly than for the transverse modes. For all GW modes from 
directions near 9 — tt, the redshift increases monotonically 
up to some limiting frequency beyond which the Taylor series 
expansion of the pulsar term which leads to Eqs. ( p7] l and ( |28| ) 
can no longer be performed. 

We now discuss the implications of this effect on the over- 
lap reduction functions. 



Relative to the longitudinal mode the redshift of vector modes 
is smaller by a factor of S and vanishes as 5 — ^ 0, but it is still 
larger than the transverse modes by a factor of 1/(5. 

The same behavior is not present in other sky locations. 
If the GW propagates perpendicular to the pulsar-Earth line 
(9 = Tr/2 + (5), then up to second order in 5 the redshifts 

^' " 2{l- S) (e"'"'^"^^'""^ - l) (longitudinal)(32) 



V. OVERLAP REDUCTION FUNCTIONS 



(1-<S) I' ') 



(shear) (34) 



are obtained. In this case for small 8 the exponential cannot be 
expanded unless /L ^ 1. For this sky location the redshift is 
always an oscillatory function of ]L. The pulse comes across 
different phases of the GW as it propagates toward Earth. 



As discussed in Section III the overlap reduction function 
for the two pulsars in Fig. [Ijis equal to 

r(|/|) = — V / df^(e2"^^^(^+^-'^^^ - 1) (35) 

= r+(|/|) + rx(|/|) + rb(|/|) + r,(|/|) (36) 
+ r.(|/|) + r,(|/|) 

where all possible GW polarizations are allowed. It is ad- 
vantageous to consider each term in the sum ([35j separately 
since various gravity theories may have different polariza- 
tion content |3 - 13, 20 1. The overlap reduction function has a 
closed analytic form for transverse GWs. The overlap reduc- 
tion function for the plus mode has been calculated by I27l 
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FIG. 5: Hellings and Downs fTl] first showed that for 
general relativity, an isotropic stochastic background of GWs 
is expected to produce the correlation shown in blue. The 
correlation for the transverse breathing mode appears in 
black. 



and is given by 

r+(e) = 3 



1 — cos ^ 



log 



1 — COS ^ 



,(37) 



where ^ is the angular separation of the pulsars. For the scalar- 
breathing mode, a closed form is given by 11211 : 



rb(0 = ^(3 + cose) 



(38) 



For the case of non-transverse GWs, the overlap reduction 
functions cannot be integrated analytically and we calculate 
them numerically. 

In general relativity the pulsar term can be excluded from 
the integral ( |35| l without any significant loss of optimal- 
ity ||24| . The reason for this is that the smallest frequencies 
that PTAs are sensitive to are ~ 0.1 yr^^, and the closest PTA 
pulsar distances are ~ 100 ly, so that fL > 10. This is shown 
in Fig.|6] where we plot the overlap reduction functions T{fL) 
with (solid curves) and without (horizontal dashed lines) the 
pulsar term for several pulsar separation angles ^ and GW po- 
larization modes. The frequencies that PTAs are sensitive to 
are to the right of the vertical dashed line at fL = 10 in each 
plot. As seen in Fig. ^3), r+(/L) is roughly independent 
of frequency over the range of frequencies relevant to pulsar 
timing experiments. The same is true for the scalar-breathing 
mode, which is shown in Fig. |6|b). It is worth pointing out 
that both r_|_(/i) and Ti,{fL) are normalized to unity for co- 
aligned pulsars. Note that the overlap reduction functions for 
all other modes are normalized with the same factor of 3/47r 
used in the +-mode. 

In Fig.|6jc), we plot the overlap reduction function Ty{JL) 
for the vector-y mode. Over the range of relevant frequen- 
cies, Ty{fL) is frequency independent for most of the pulsar 
separation angles shown. For co-aligned pulsars, however, 
Ty{fL) retains frequency dependence well into the range 
of pulsar timing frequencies, and takes on values an order 
of magnitude higher than those obtained by T^{fL) and 
TbifL). 



Similar behavior is shown in Fig.|6|d), where we have plot- 
ted the overlap reduction function for the scalar-longitudinal 
mode. Here Ti{fL) retains frequency dependence through- 
out the relevant frequency range for each of the pulsar separa- 
tion angles shown. For the case of co-aligned pulsars, Ti{fL) 
does not converge, and for separation angles that do converge 
F; {fL) takes on values that are at least an order of magnitude 
larger than those obtained by r+(/L) and Tb{fL). 

For co-located pulsars we can understand the behavior of 
the longitudinal mode analytically. In the problematic sky 
region {6 w tt), Ti{fL) is proportional to the square of the 
redshift. 



TiifL) cc 2tt 



cos 



,-27ri/L(l+cos( 



rf(cos ( 



1 



(39) 



4(1 + cos 61)2 

which may be evaluated analytically. In the limit of large fL, 

TiifL) = 7r{37/6-47 -l/(7r(/L)2)+4Ci(47r/L) 

- 4 log (47r/L) + Si(47r/L)} (fL > 1) 

- (37/6 - 47) TT - 47r log (47r/L) + n^fL, (40) 

where 7 is Euler's constant. The overlap reduction function 
Vi{fL) is roughly proportional to fL in this limit. Eq. ( |40l ) 
is shown along with the numerically integrated overlap reduc- 
tion functions in Fig. |6|d) and, with the exception of the sin- 
gular behavior near the origin (where the large fL approxi- 
mation is not valid), agrees well with the numerical Ti{fL) 
curve for co-ahgned pulsai^s (^ — 0). 



VI. OVERLAP REDUCTION FUNCTIONS FOR THE 
NANOGRAV PULSARS 

The NANOGrav PTA consists of 24 pulsars. The Australia 
Telescope National Facility (ATNF) data for the distances to 
these pulsars is given in Table |l] [2_8J. Using a simple nu- 
merical integration scheme, the overlap reduction function for 
each pulsar pair was computed. The main difference relative 
to the previous section is that we are including the effect of 
different pulsar distances. Results are given in Fig. |7] (a)- 
(d) and show that the calculated values of r(/) are consis- 
tent with the more simple results discussed in Section [V] for 
the non-transverse modes for frequencies up to ^ 10~ Hz. 
Pulsar pairs with the smallest (^ < 12°) separation angles 
(starred curves in Fig.|7](b), (d)) for non-transverse polariza- 
tion modes are characterized by large values of the overlap 
reduction function and monotonic growth up to some limit- 
ing frequency. Pulsar pairs with larger (^ > 12°) separation 
angles (un-starred curves in Fig. |7] (b), (d) and all curves in 
Fig.|7]) do not display monotonic growth up to a limiting fre- 
quency, but still result in much larger values than those of the 
plus and cross modes. Fig. [T] shows that sensitivity is greater 
for scalar-longitudinal and vector modes than for the tensor 
and scalar-breathing modes, and increases rapidly for pulsars 
that are nearly co-aligned in the sky. Theoretical sensitivity 
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FIG. 6: (color online) T{fL) with (solid curves) and without (horizontal dashed lines) the pulsar term for the various 
polarization modes: plus (a), breathing (b), shear (c) and longitudinal (d). In the latter two modes, smaller pulsar separation 
angles are characterized by retained frequency dependence in T{fL) in the range of frequencies relevant to pulsar timing 
experiments. Nearly all the non-transverse curves eventually converge, but at rather high values of T{fL) relative to the 
transverse modes, indicating increased sensitivity to GWs with these polarizations. We have plotted the large limit 
approximation (|40]i (dashed black curve) along with Vi{fL) in (d), which is in good agreement with the ^ = curve. 
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TABLE 1: NANOGrav Pulsai- Data 



estimates using (|40J support this result (see Appendix |B). In 
these plots, increased variation in curve amplitudes appears 
due to the fact that the pulsars are not equidistant from the 
Earth. 

Over the entire range of frequencies plotted for pulsar tim- 
ing experiments (between ^ 10^^ and ^ 10^^ Hz), the over- 
lap reduction functions are approximately constant. In prac- 
tice, some optimality will be lost due to the fact that pulsar 
distances are known at best to only ^ 10% | 



VII. DISCUSSION 

Direct detection of GWs might be possible in the next 
decade using a pulsar timing array. A detection would pro- 
vide a mechanism for testing various metric theories of grav- 
ity. To develop optimal detection strategies for stochastic 
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FIG. 7: r(/) for some of the NANOGrav pulsar pairs. Pulsar pairs, along with their angular separation in degrees, are shown 
with each curve. As / increases, r(/) approaches a constant value. The asterisk indicates the NANOGrav pulsar pair with the 
smallest angular separation 3.35 degrees). Note the larger values of the r(/)s for this pair. 



backgrounds in alternative theories of gravity, we have com- 
puted overlap reduction functions for all six GW polarization 
modes, including four modes not present in general relativity. 

We began by introducing the redshift induced by GWs 
of various polarizations, along with the polarization tensors 
unique to each mode. We then derived the optimal detection 
statistic for an unpolarized, isotropic stochastic background 
of GWs, and defined the overlap reduction function as a ge- 
ometrical quantity that is proportional to the expected cross 
correlations. 

We examined the redshifts induced by GWs of various po- 
larizations on the pulsar-Earth system. We show that when 
the GWs are coming from roughly the same direction as the 
pulses from the pulsar, the induced redshift for any GW polar- 
ization mode is proportional to fL, the product of the GW fre- 
quency and the distance to the pulsar. When the GWs and the 
pulse direction are exactly parallel the redshift for the trans- 
verse and vector modes vanishes, but it is proportional to fL 
for the scalar-longitudinal mode. The vanishing contributions 
from the tensor, vector and scalar-breathing modes are there- 
fore not a result of the pulse surfing the GW. Rather, it is the 



transverse nature these modes that is responsible for the van- 
ishing response. In this limit we also show that the redshift is 
proportional to the relative velocity of the pulsar-Earth system 
(Lh), which is the same when the pulse is emitted and when 
it is received. 

We find that the overlap reduction functions for non- 
transverse GWs are characterized by frequency dependence 
that is significant for nearby pulsar pairs. The values of the 
overlap reduction function increase by up to one order of mag- 
nitude for the vector polarization modes and up to two orders 
of magnitude for the scalar-longitudinal mode. Pulsar timing 
arrays are significantly more sensitive to scalar-longitudinal 
and vector GW stochastic backgrounds. 

Next, we used current pulsar distance and sky-location data 
from the ATNF pulsar catalog to calculate the overlap reduc- 
tion functions for each pulsar pair in the NANOGrav pulsar 
timing array. Over the range of frequencies relevant to pulsar 
timing array experiments, these overlap reduction functions 
for all polarization modes are roughly constant for most pul- 
sar pairs. For nearly co-aligned pulsars, the overlap reduc- 
tion functions for scalar-longitudinal and vector modes ex- 
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hibit marked frequency dependence and asymptote to much 
larger values than the overlap reduction functions for trans- 
verse modes. In fact for a pair separated but about 3° we find 
a sensitivity increase of about a factor of 10^ for longitudinal 
modes. 

The results discussed here may be compared to other re- 
cent work. Lee et al. |21| calculated the cross-correlation 
functions for stochastic GW backgrounds including all six 
GW polarizations, and found that the correlation functions 
for non-transverse GWs are frequency dependent, as well 
as an increased response in the cross-correlation to scalar- 
longitudinal GWs, in agreement with our results. This 
work was done in the context of the coherence statistic fTP\ 
for stochastic background detection, rather than the optimal 
statistic 1241 • The coherence statistic is a measure of good- 
ness of fit of the pulsar-pair cross-correlations to the Hellings- 
Downs curve. For non-transverse modes there is no Hellings- 
Downs curve because the overlap reduction functions remain 
frequency dependent for large fL. Lee et al. solved this 
problem by simulating GW backgrounds and finding effec- 
tive background-dependent Hellings-Downs curves for these 
theories. In the context of the optimal statistic this is a non- 
issue: The frequency dependent overlap reduction functions 
can be used to construct the optimal filter in Eq. (|22]). This 
is identical to what is done for LIGO stochastic background 
optimal filter construction |25|, where the overlap reduction 
functions are also frequency dependent. 

Alves and Tinto ll22l have estimated antenna sensitivities 
to GWs of all six polarization modes by assuming a signal- 
to-noise ratio of 1 over 10 years and calculating the noise 
spectrum. Their results indicate an increase of two to three 
orders of magnitude in sensitivity to scalar-longitudinal mode 
GWs compared to that of plus and cross mode GWs. To ex- 
plain this effect Alves and Tinto compare the effect of a tensor 
GW propagating orthogonally to the pulsar-Earth system, and 
a scalar-longitudinal GW propagating in a direction parallel to 
the pulse direction. They argue that the increased sensitivity to 
longitudinal GWs is due to the amount of time a longitudinal 
GW affects the pulsar-Earth radio link. 

We have compared the effect of GW propagation from di- 
rections near the pulsar and orthogonal to the pulsar-Earth 
system for all polarization modes. For GW propagation di- 
rections parallel to the pulse direction we find that the red- 
shift induced by a gravitational wave is large, and seemingly 
divergent when the GW and pulse directions are exactly par- 
allel. This apparent divergence occurs for longitudinal, trans- 
verse, and shear modes alike. In that limit, however, the diver- 
gent term in the redshift that comes from the relationship be- 
tween time and affine parameter derivatives cancels because 
the phase of the GW pulse when pulse is emitted is nearly 
equal to the phase of the GW when the pulse is received (see 
Eqs. ( |A9| l, ( |23] l and (p4]i). The redshift becomes proportional 
to the relative velocity of the pulsar-Earth system and a mode- 
dependent geometrical projection factor for all GW polariza- 
tion modes. In this limit the relative velocity of the pulsar- 
Earth system is approximately equal when the pulse is emitted 
and received. For transverse and shear modes the projection 
factor vanishes when the GW and pulse directions become 



parallel. For longitudinal modes the geometrical factor goes 
to a constant, so that the pulsar-Earth system is very sensitive 
to GWs from directions near the pulsar This is the physical 
origin of the increased sensitivity to scalar-longitudinal GWs. 
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Appendix A: Analog to Detweiler's equation for vector and 
scalar polarization modes 

Using the techniques used previously by 1*23^, '241 for Ein- 
steinian GWs, we compute the induced frequency change for 
non-Einsteinian GWs. We begin by considering the metric 
due to a longitudinal mode gravitational wave perturbation: 



gab = Vab + habit - z) 

f ~l \ 
10 
1 

V 1 + hL J 



(Al) 



Given a null vector = — a, — /3, — 7) in Minkowski 
space (where a, f3, 7 are directional cosines) the correspond- 
ing perturbed null vector is given by 



/ 1 

—a 

v-7(i-¥)y 



(A2) 



From the geodesic equation, the t-component of cr" must sat- 
isfy 



da* 



-1 nhO^ 



(A3) 



where 



1 . 



(A4) 
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Now we may write the geodesic equation as 

1. „ . 



dcr* 



To zeroth order in hj^. 



2 

1 • 



(A5) 



hL 
2 



allowing us to write the geodesic equation as 

da* dv 1 



dX 



(A6) 



(A7) 



We now need to express the time derivative of the metric per- 
turbation, hi,, as a derivative of the affine parameter A. Since 
/i/, = hi,{t — z), we may write 



dhL 
~dX 



dh^ dt dhL dz 



dt dX 
dhL dt 



dz dX 
dh^ dz 



dt dX dt dX' 



(A8) 



Identifying the relations ~ ^ ^^'^ 3X ~ we obtain 

the relation 



hr^ 



dhi 



1 



dhj 



dt v{l + i) dX 
which makes the geodesic equation 



dv 
dX 



dhr 



1^- 2 2 1 t^7^ 



Integrating both sides, we obtain 



exp 



1 1' 
2(1 + 7) 



Ah, 



(A9) 



(AlO) 



(All) 



where A/i; = hf — h^. Expanding to first order in /i^, we 
may write 



7 



1 

"2(1+7) 
^„o2 a 



Ahr 



2(1 + COS 6'p) 



Ahr 



(A12) 
(A13) 



The derivation for vector modes is nearly identical to that of 
the longitudinal mode. For the sake of brevity we only detail 
the vector-y mode in the remainder of this document. For the 
vector-y mode, the metric perturbation takes the form 



gab 



f -1 \ 
0100 
I hy 

V hy 1 J 



(A14) 



The null vector becomes 



V 



¥-7 / 



(A15) 



Following the same algebraic steps used above, one obtains 
the geodesic equation 



which leads to 



da' 
dA 

dv 
dX 



dv 
dX 

v"fP dhy 

"(1 + 7) dT' 



(AI6) 



(A17) 



Integrating this expression and expanding the result to first 
order in Ahy produces the result 



cos Up sm <pp sill t'p 



(1 + cos 9p) 



(A18) 
(A19) 



where Ahy — hy — hP. 

For comparison, we also include the results for the plus, 
cross, vector-x, and breathing modes. For the plus mode, we 
obtain 



2(1+7) 



Ah. 



(A20) 



, sm Op + sm Op sm 



2(1 + cos 6'p) 



-A/14 



for the cross mode, 

Ve - Vp 



af3 



for the vector-x mode, 

Ve - Vp 



(1+7) 



Ah, 



cos 0„ sm t/„ sm < 



(1 + cosOp) 



-Ah, 



a7 



-Ahr 



v„ 



(1 + 7) 

cos 0„ cos (j)„ sin t 



(1 + cos 6'p) 



-AK 



and for the breathing mode, 

Ve - Vp {a^ + 



Ahh 



2(1+7) 

ip sill wp 



cos^ (hn sin^ 



sin^ Or, sin^ 1 



2(1 + cos 6'p) 



(A21) 

(A22) 
:; (A23) 

(A24) 
; (A25) 

(A26) 

'^Ah,. 
(A27) 



Here, Ahj^ = h^ — h^, and we can identify these expressions 
with (|2]i. 
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Appendix B: Sensitivity bounds for Ti{f) 

In this section we derive an estimate for the value of r;(/) 
for pulsars that are nearly co-aligned in the south-pole region 



of the sky. Consider a two-pulsar Earth system in which the 
pulsars are separated by angle ^ and are equidistant from the 
Earth (Li — L2). In this scenario the overlap reduction func- 
tion for a longitudinal mode GW is given by 



T if) f dO, '"''^^'^ ^ (cos^cos^ + sin^sin^cos^)^ 
' Jg2 2(1 + 0086* (1 + cosScos^ + sin0sin^cos( 



\ I g — 27ri/L(l+cos COS ^+siii ^ sin ^ COS 0) ^ 



r 



(Bl) 



We will make the following approximations: suppose that 
6 = TT — S, such that S < 1 (this is necessary to achieve 
a large enough response in cross-correlation). Doing the rele- 
vant series expansion in S and retaining terms to second order, 
the integrand simplifies significantly: 



Identifying S ^ where n is a real constant satisfying n E 
(0, 1], we may now re-express r;(/) as a function of ^: 



Tiif) = [nfLf / # / d5(l-5^) 

JQ Jo 



X (l - (5^ - (5i^cos( 



(B2) 



Here we have changed the bounds of integration for 6 based 
on the fact that the argument of the exponential terms in ( |B1| | 
is necessarily small. 

Completing the integration over (p eliminates factors of ^ 
in the remaining integrand and the remaining mathematics is 
trivial: 

Tiif) = 27r(^/i)2 f'^d5{l~5^f 
Jq 



2n 



(7r/L)3/2_-y^ 



(B3) 



r/(/) 



27r 



2 1 



«)3 3 «) 



(B4) 



Some potential values of Ti{f ) are given in Table [n] For ex- 
ample, for a pulsar pair that is separated by only ^ 3°, with 
ri = 1 we obtain F; (/) ^ 10"*, which is of the same order of 
magnitude as that obtained for the J1853-J1857 pulsar pair in 
Fig.|7];b). 
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